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Abstract 

Stationary bifurcations in several nonlinear models of fluid conveying pipes fixed at 
both ends are analyzed with the use of Lyapunov- Schmidt reduction and singular- 
ity theory. Influence of gravitational force, curvature and vertical elastic support on 
various properties of bifurcating solutions are investigated. In particular the condi- 
tions for occurrence of supercritical and subcritical bifurcations are presented for 
the models of Holmes, Thurman and Mote, and Paidoussis. 
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Introduction 



The nonlinear dynamics behavior of pipes conveying fluid displays in- 
teresting and paradigmatic behavior with important practical implications and 
it offers a rich setting for development and testing of nonlinear dynamics the- 
ory. Bifurcation theory represents one of the main subject areas of nonlinear 
dynamics, analyzing behavior either in the vicinity of trivial solutions (local 
bifurcations), or the existence in-the-large of connected sets of non-trivial so- 
lutions (global bifurcations). In the study of local bifurcations, particularly 
interesting is the influence of various parameters figuring in the governing 
equations on the location and stability of fixed points and on classification of 
bifurcations as supercritical or subcritical. In a supercritical bifurcation there 
is no discontinuous change in size and form of the attractor, and after the bi- 
furcation the new (enlarged) attractor contains within itself the old attractor. 
On the other hand, in subcritical bifurcation the attractor disappears followed 
by a jump of the system to a remote and completely new attractor, via a fast 
dynamic transient. For the complete understanding of the dynamics, both 
stationary and nonstationary (dynamic) aspects of the bifurcation theory are 
important. 
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In the case of fluid-conveying pipes, nonstationary analysis is related to 
the study of supercritical and subcritical Hopf bifurcations, the determination 
of the amplitude associated with flutter and the dependence of oscillation fre- 
quency on the amplitude. Stationary analysis, on the other hand, is concerned 
with stationary bifurcations i.e. changes in the equilibrium point structure 
of the underlying equations, which due to reflection symmetry are of pitch- 
fork type. Furthermore interplay between stationary and dynamic bifurcations 
may lead to more complicated dynamics such as spatio-temporal intermittency 
and chaos (Argentina and Collet, 1998). In spite of many publications devoted 
to these subjects, a thorough and systematic analysis of nonlinear models of 
fluid conveying pipes, particularly of pipes supported at both ends, and re- 
lated bifurcations is still lacking. Following the pioneering works by Holmes 
(1977, 1978) further work of interest appeared in Ch'ng (1977, 1979) and Lunn 
(1982). As a paradigm of the studies performed so far, the work of Holmes 
based on the center manifold analysis concentrated on one-equation model 
which does not take into account any gravitational nor tensioning effects. The 
coefficients figuring in the normal form of the obtained bifurcation equation 
were given only as numerical values, disguising the influence of involved phys- 
ical quantities. Finally, the complete two-equation model of Paidoussis (2003) 
has only recently been introduced in its correct form and it has never been 
used. A wealth of up-to-date information related to the published works on 
nonlinear aspects of fluid-conveying pipes with supported ends may be found 
in Paidoussis (1998, 2003). 

In order to improve upon the current status and in order to prepare 
the ground for the study of spatio-temporal intermittency (and chaos), we 
study several well-known and frequently used nonlinear models in which we 
focus on nondegenerate local bifurcations from the trivial solution to another 
stationary nontrivial solution. In particular, the object of our analysis are pipes 
fixed at both ends shown in Fig. 1 and nonlinear models of Holmes (1977), 
Thurman and Mote (1969) and the complete nonlinear model of Paidoussis 
(2002), [the latter being based on the model of Semler at al. (1994)]. The 
manner of presentation is such that each model is considered as a special 
case of the complete nonlinear model, revealing the influence on the complete 
dynamics of different linear and nonlinear terms and quantities. Moreover the 
increasing complexity of each model leads to a better understanding of the 
complete nonlinear model. Each model is presented in a separate section of the 
paper, with subsections related to certain important parameters influencing 
the dynamics, such as the gravitational force or pipe curvature. 

Since this exposition relies strongly on the Lyapunov-Schmidt reduc- 
tion, a brief description of the procedure is presented in Appendix A. The 
interested reader may be find extensive treatment in books (Golubitsky et 
al., 1985; Golubitsky and Schaeffer, 1988). An important characteristic of the 
Lyapunov-Schmidt reduction is that the procedure for obtaining normal form 
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of the bifurcation equation is pursued by analytical means, without using any 
numerical or truncating procedures. Hence the normal form is exact and influ- 
ence of each parameter on the type of bifurcation may be traced and analyzed. 
An important insight gained by using this method is reflected in the fact that 
we obtain exact analytical solutions in the vicinity of bifurcation point for each 
model and derive exact conditions that classify bifurcations as supercritical or 
subcritical. 



1 The Model of P. J. Holmes 



1.1 Bifurcations in the stationary model without gravitational effects 



Holmes considered pipes with supported, non-sliding ends and obtained 
a nonlinear equation of motion by adding to a linear equation the nonlinear 
term corresponding to the mean, deformation- induced tensioning. The com- 
plete equation of motion in nondimensional form reads 



tttt 

av 



+ v "" - (r - p 2 - X -A \v'f + aA (v 1 | v})v" (1) 

-7([1 -S]v'y + 2y/ppb' + ai)+ v = 0, 

where v = v(£,t) denotes the lateral deflection normalized by the length of 
the pipe. In addition, T represents the tensile force on the pipe, (3 represents 
the mass ratio, p the flow velocity, A the axial stiffness, a is related to the 
viscoelastic structural damping, a represents fluid damping, and 7 denotes 
gravitational effects. Explicitly, 

/ EI \ 1/2 a M M + m T3 

a= {wT^) V> ^ = MW 1 = ^T L9 > (2) 

where a is the coefficient of Kelvin- Voigt damping in the pipe material, M 
and m are the mass of fluid and pipe per unit length, respectively, and L is 
the length of the pipe, while 

r _ T L 2 _ EAL 2 cL 2 

EI' EI ' [EI(M + m)}V 2 ' { ' 

In the above equations T , EI, L, A, I, E and c are the longitudinal externally 
applied tension, flexural rigidity of the pipe, length of the pipe, cross-sectional 
area, area-moment of inertia, modulus of elasticity and constant of damping 
(due to friction), respectively. Moreover v' = dv/d£, v = dv/dt. With £ = s/L, 
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\v I denotes the L 2 norm 

W\ = ([ 1 v> 2 (Od0 1/ \ (4) 
Jo 

while the (v' \ v ^ denotes the L 2 inner product 

(v , \v) = £v , (t)v(t)d£. (5) 

The boundary conditions are v — v" — at £ = 0, 1. The corresponding 
stationary equation, neglecting gravitational effects, applied tension and non- 
linear dissipative term following a Lyapunov-Schmidt reduction, locally may 
be put into one-to-one correspondence with solutions of a single algebraic 
equation 

g(p, p) = (J(mr) Vl) /i 3 - (utt) 3 pp + 0(3) = 0, (6) 

indicating that the corresponding bifurcation is always a supercritical pitch- 
fork one, since 

hnir) A A> 0. 
4 



Using a two-mode Galerkin discretization applied to the complete equa- 
tion (1) and assuming a — g — V — 0, Holmes (1977) obtained numerical val- 
ues for terms multiplying Ap? and pp respectively in the normal form equation 
(6). We have shown here that the sign of these terms is unalterable and we 
have presented their analytical form. Applying the Lyapunov-Schmidt reduc- 
tion to the same stationary equation and taking applied tension into account, 
the following normal form of the bifurcation equation is obtained: 



g(p,p) = - 4 (nn + £^yA 



sin' 



(-) 

\ nn ) 



^ 3 -(mr) 3 pp + #(3) =0. (7) 



In this case too, the type of pitchfork bifurcation is supercritical and may not 
be altered, irrespective of the value of tension T or axial flexibility A. 



Since the influence of gravitational effects was not considered in the 
original treatments (Holmes, 1977, 1978), nor in subsequent studies, we ana- 
lyze this effect along with the effect of tensioning in the next section, again 
restricting analysis to the stationary version of equation (1). 



4 



1.2 The role of gravitational effects 



The stationary equation representing nonlinear dynamic behavior of 
a pipe conveying fluid including gravitational and tensioning effects may be 
written in the form 

v"" -{p 2 -T- l -A \vf)v" - 7 ([1 - £]„')' = 0, (8) 

where v, as in (1) denotes the displacement in the z-direction. The boundary 
conditions are v = v" = at £ = and £ = 1. Introducing the variable 
a 2 = p 2 — (r + 7), integrating once and setting dv/d^ = u we obtain the 
equation 

u" + (a 2 +rt)u-^A\uf = 0, (9) 

with boundary conditions tt'(£ = 0) = = 1) = 0. The one-dimensional 
kernel of the corresponding linear operator 

L ■= u" + (a 2 + 7f )u, 

is now spanned not by a sine function, as in the case without gravitational 
effects, but by Bessel functions. The proof of this proposition requires trans- 
formation of the equation 

u" + (a 2 + 7 £)m = 0, (10) 

to the equation 

t 2 w" + tw' + (t 2 -(-) 2 )w = 0, (11) 
3 

where in the above equation prime denotes differentiation with respect to t. 
First, substitution (a 2 + 7^)7^ = y, o = —2/3, leads to the equation 

u" + yu = 0, 

where the prime denotes differentiation with respect to y. The subsequent 
substitution t = |?/ 3 / 2 , u = ^Jyw leads to equation (11). This equation may 
be solved in terms of Bessel functions so that the corresponding solution of 
the linearized eq. (8) is 



v(0 = (a 2 + 7 



J 2 /3 [~ 3 -) J-2/3 ^3 

2a 3 \ (2{a 2 + 



-J-2/3 g~ J2 / 3 



3 7 ) \3 7 / 



(12) 



For gravity and fluid motion in the same direction, and assuming T = 0, 
the location of bifurcation points (nontrivial equilibria) is obtained from the 
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expression 

p n = nn + -^- + 0( 7 2 ), n = 0,l,2,.... (13) 
4 nn 

Details of the derivation of equations (10) and (13) are presented in Appendix 
B. The above expression suggests that bifurcation points are shifted as com- 
pared to the case without gravity in the direction of increasing velocity by 
the amount 7/4n7r. Correspondingly, bifurcations occur for higher velocities 
compared to the case when gravity is not taken into account. As the velocity 
increases, the effects of gravity are weaker, so that finally locations of bifurca- 
tion points are the same as for the case without gravity. This is an intuitively 
pleasing result as one expects high enough velocity to annul the effects of 
gravity. If T is not neglected, the expression for critical velocity is 

Pn = mx + -^- + \— + 0(2), (14) 
4 nn 2nn 

= n7T + Il±^E + tf(2), n = l,2,... 
4 nn 

The resulting bifurcation diagrams are presented in Fig. 2. Hence the roles 
of tensioning and gravity are similar; however, as the velocity increases the 
influence of gravity diminishes twice as fast as tensioning effects. The other as- 
pect of this phenomenon may be observed by considering the distance between 
bifurcation points, in units of velocity, given by the following expression: 

i 7 i r 

When 7 = T = 0, bifurcation points are equidistant from each other, this 
distance being equal to n. When 7 and T are not equal to zero, the distance 
is an increasing sequence, whose limit {n — > 00) is equal to n. 



If the direction of the fluid is opposite to the direction of gravity and 
assuming no tensioning effects, bifurcation points are located at: 

p n = nn-\^- + 0( 7 2 ), n = 0,l,2,.... (16) 
4 nn 

In this case bifurcations occur for velocities smaller than in the case without 
gravity, and again for high enough velocities the effects of gravity on the 
location of bifurcation points are annulled. Inclusion of tensioning effects yields 

1 7 1 r 1 7 - 2r 

Pn = ™r - + - — + 2 = nn - -± + 0(2), n = 1, 2, 

Ann 2nn 4 nn 

and the corresponding bifurcation diagrams are presented in Fig. 3. The 
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matching bifurcation interpoint distances are 




i 7 i r 

Ann(n + 1) 2ir n(n + 1)' 



In this case there is a competition between gravity and tensioning effects, with 
the possibility that they annul each other when 7 = 2I\ 

Equation (11) yields an approximate solution 



which for 7—^0 approaches the solution sin(n7r£), the solution of the linear 
equation corresponding to (8) without gravitational effects. In order to shed 
more light on the physical aspect of the above results, it is instructive to call 
upon an important relationship between the solution p(p) of the algebraic 
equation g(p, p) = and the solution of the full problem f(£). Near the bi- 
furcation point the nontrivial solutions of (8) have the spatial structure of the 
basis vector u\ that spans the one dimensional kernel of the corresponding 
linear operator. Hence, the solutions /i of the reduced bifurcation equation 
g(p, p) = are related to solutions of Eq.(8) as 



so that the complete solution of (8) in the vicinity of the bifurcation point 
may be written as 



satisfying boundary conditions v (0) = v(l) = 0; v"(0) = v"{l) = 0. 

The analysis of the Holmes' model reveals several important features. 
Irrespective of the inclusion or exclusion of gravity and/or tensioning effects, 
the stationary bifurcation is of supercritical type. Although gravity and ten- 
sion have a similar effect on the location of bifurcation points, only gravity 
changes the solution locally. In the vicinity of the bifurcation point, the solu- 
tion (eigenfunction) is a sine function for 7 = and a Bessel function for 7 7^ 0. 
Moreover, both gravity and tension increase the distance between bifurcation 
points along the parameter space and this effect is particularly noticeable for 
low velocities, while high velocities compensate the influence of gravity and 
tension as intuitively expected. 




(17) 



v = pui + $(/x 2 ), 




(18) 
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2 The Model of Thurman and Mote 



2.1 Bifurcations in the stationary model 



The model of Thurman and Mote (Thurman and Mote, 1969; Paidous- 
sis, 1998), which considers both lateral and axial deflections, was derived un- 
der the following assumptions: (i) no gravity force (ii) steady flow velocity (iii) 
linear moment-curvature relationship and (iv) a simple approximation of the 
fluid velocity. For the study of the normal form of the bifurcation equations 
it is more instructive to consider this model as a special case of the complete 
nonlinear model of Paidoussis (2003) 1 , to be analyzed in the next section. 
The nondimensional equations for spatial motion of the complete nonlinear 
model are 



w"" + (p 2 - (r - u))w" + 7 w' + (r - A - n) («,v + u'w" + ^w' 2 w") 

-(3u"'w" + Au"w"' + 2u'w"" + w'u"" + 2w" 3 + 2w /2 w"" + 8w'w"w'") 



-7 



w'u' + ho' 3 -{1-0 (-w" + u"w' + u'w" + ^w' 2 w" 



0. 



(19) 



(p 2 - A)u" - (w"w"' + w'w"") - 7 \^w' 2 - (1 - i)w'w'' 

+ (r - A - n) w'w" = o. (20) 

where the dimensionless tension and the pressure at the downstream end are 

T(L)L 2 n _ P(L)L 2 



r = 



EI 



EI 



21) 



respectively, while u denotes the dimensionless longitudinal deflection in the x 
direction (direction of gravity), and the other quantities being the same as in 
the model of Holmes, with w here replacing v. A prime denotes differentiation 
with respect to the Lagrangian variable x which may be used interchangeably 
with £. The following assumptions may be attributed to the model of Thurman 
and Mote: 



7 = 0, 



n = o 



T L 2 
EI ' 



K 



o, 



where T =const, denotes externally applied tension and n 2 denotes curvature. 
Consequently the equations of motion of this model are 



The corrected version of equations is given in Appendix T.4. 
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w »» + ( p 2 _ r)w // + (r-A) {w"u' + uV + |u/V') = 0, 

(p 2 -^K + (r-^) w V = o. (22) 

The boundary conditions are u(0) = w(0) = u(£ = 1) = w(£ = 1) = 0, with 
the additional condition w"(0) = u>"(£ = 1) = for a simply-supported pipe 
or w'(0) = = 1) = for clamped-clamped one. 



Fixed points, as in the model of Holmes for 7 = 0, are located at 
positions 

Pn = ™ + \ — + tf(r 2 ), n = 0,l,2,.... 
I nix 

as easily seen from the following argument. The kernel of the linear operator 
L corresponding to the system consisting of Eqs. (19) and (20) is obtained by 
solving the equation 




<' + (p 2 - r)< 



V 



Ur 



(23) 



where 



w {x) 



\U {X) 

satisfying boundary conditions 



Eoo 
n=l 



sinn7r£, 



(24) 





Hence, the position of bifurcation points is readily obtained from equations 
(23) and (24). A straightforward calculation indicates that the kernel of L is 
spanned by 



KerL( 



(w \ 


)=R< 






■ e \ 

sin nnt ) 


,p 


! 




\u ) 






w 


i 



(25) 



Calculation of terms defined in the Appendix A, leads to the following bifur- 
cation equation 



g(p,p) = -(n7t) 4 [a({3-3)} /i 3 - {rmf pp + $(x 3 ) = 0, (26) 



where 



a = r-A, 



r-A 



p 2 -A' 

In the vicinity of the bifurcating solutions, the solutions of the complete model 
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have the form 

" iC) =MP) 1 | sin for* • (27) 

where p(p) represents the solution of equation (26). Hence, again the bifurca- 
tion is of the pitchfork type; however, there are important differences between 
the model of Holmes and the model of Thurman and Mote. It is immediately 
evident that bifurcation may be of supercritical pitchfork type for 

a(/3-3) > 0, 

and of subcritical type, shown in Fig. 4, for 

a((3-3) < 0. 

Writing explicitly the expression on the left side of these inequalities as 

(r - A) 2 



3(r-^i), 



and recalling that p 2 = T + (mr) 2 , implies that two cases may be analyzed, 
depending on whether \a\j p 2 n < 1 or p 2 n j\o\ < 1. However, since only the 
first point of instability corresponding to n — 1 is relevant from the physical 
point of view the former condition is of no practical importance although it is 
important for understanding the mathematical aspects of the model. 

Case 1: \a\j p 2 n < 1. This situation corresponds to high fluid veloci- 
ties (large n). In this straightforward calculation yields the following 
conditions for the occurrence of supercritical and subcritical bifurcations: 



Supercritical condition: 



T < A. (28) 



Subcritical condition: 



T > A. (29) 
Based on (3) the above inequalities may be also expressed as 

T < EA, 

and 

T > EA. 

Case 2: p 2 n l\oi\ < 1. This case is applicable to low fluid velocities (small 
n ) and large |T — A\ , so that either a highly flexible pipe is considered or the 
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effects of tensioning are large. Conditions for the occurrence of supercritical 
and subcritical bifurcations are 



Supercritical case: 




(30) 



Subcritical case: 




(31) 



For low fluid velocities, both conditions strongly depend on values of T and A, 
since n is small so that supercritical bifurcation is more likely in a pipe of high 
axial flexibility and low externally applied tension while subcritical bifurcation 
is more likely in a pipe with small axial flexibility and high externally applied 
tension. 

Hence, an important feature of the Thurman and Mote model is depen- 
dence of the bifurcation type (supercritical or subcritical) on the relationship 
between elastic characteristics of the pipe and externally applied tension, with 
the velocity-dependent term figuring only in inequalities relevant at low fluid 
velocities (the case of physical validity). At this point it should be mentioned 
that if p 2 does not appear (Eq. (5.62) of Paidoussis (1998)) in the equation 
for u, only one classification condition is obtained in the form of inequalities 
(28) and (29); thus from the point of view of the bifurcation theory this ver- 
sion of the model assumes high fluid velocities. On the other hand, since the 
case of high fluid velocities requires large n and is therefore of no practical 
importance, the term p 2 is necessary for the model to be useful in bifurcation 
analysis. Comparison with the model of Holmes, which assumes just lateral 
deflections and for which only supercritical bifurcation is possible, reveals that 
inclusion of the equation for axial deflections makes subcritical bifurcation also 
possible. 

2.2 Thurman and Mote model with curvature 

The inclusion of the curvature term EIk 2 in the model of Thurman and 
Mote may again be considered as a special, reduced version of the complete 
nonlinear model. Retaining simplifications of the Thurman and Mote model 
in effect, namely 7 = 0, n = 0, Y = T L 2 /EI, k 2 = 0, T =const., but 
keeping terms that arise from the curvature effect, the following dimensionless 
equations of motion are obtained: 



-(3u"V + Au"w"' + 2u'w"" + w'u"" + 2w" 3 + 2w' 2 w"" + 8w'w"w'") = 0. 



w »» + ( p 2 _ r)u /' + (r-A) {w"u' + u'w" + -w' 2 w") (32) 
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(p 2 - A)u" - (w"w"' + w'w"") + (V-A) w'w" = 0. (33) 

Boundary conditions are the same as in the basic Thurman and Mote model. 
Arguments used in previously discussed models lead to, as expected, same 
locations of fixed points 

Pn = rm + \— + tf(r 2 ), n = 0,l,2,.... 
2 mi 



Evaluation of terms of the bifurcation equation shows that it has the 
following form: 

g(x,p) = ^rmf[(7 - P){rm) 2 - ± a (f3 + 3)]x 3 - (nvr) 3 xp, (34) 

where 

a = T — A, g= £z4±W = 1+ 



p\ - A a + inn) 2 ' 

Hence the bifurcation is of supercritical type if 

(7-P)(nn) 2 -±a(f3 + 3)>0, 



and subcritical if 



(7-/3)(mr) 2 - i«(/? + 3) < 0. 



More insight into these inequalities is gained by considering specific cases 
dependent on the fluid velocity. 

Case 1: \aj p 2 \ < 1. This case corresponds to high fluid velocities (large 
n). A straightforward calculation leads to conditions for the occurrence of 
supercritical and subcritical bifurcations: 

Supercritical condition: 

r<^+^(ri7r) 2 . (35) 



Subcritical condition: 

r>^+^(n7r) 2 . (36) 

A very large and hence a nonphysical value of n would be required to make 
the term (mr) 2 dominant in the above inequalities. Since with increasing fluid 
velocity p the effective stiffness of the pipe diminishes, the model strongly 
prefers supercritical bifurcation. However, since only the first mode is relevant 
from the physical aspect, the interest is in the case below. 
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Case 2: |(mr) 2 /o:| < 1. This case may be attributed to low fluid veloc- 
ities (small n), and is therefore of physical and practical interest. The corre- 
sponding conditions are: 



Supercritical condition: 

7 



r <A+-(nir) 2 . (37) 



Subcritical condition: 

r>„4+^(wr) 2 . (38) 

The inclusion of the curvature term in the model increases the effective axial 
flexibility (the term A + |(w) 2 ) as compared to the basic Thurman and Mote 
model (the term A— ^(nir) 2 ), so a slight preference is given to the supercritical 
bifurcation. Taking into consideration that n is small, tensioning effects may 
be dominant, for example in the case of short pipes with high flexural rigidity. 



2.3 The effect of elastic support 



In order to investigate the essential features of the added elastic support 
which involves distributed springs along the length of the pipe, we use the 
least complex model, the model of Thurman and Mote model without gravity 
effects. The corresponding equations are 



w 



+ (p 2 - r>" + (r - A) {w"u' + u'w" + >'V) + Kw = 0, (39) 



(p 2 - A)u" + (T-A) w'w" = 0. (40) 

The algebraic bifurcation equation is identical to eq. (26), hence elastic support 
does not change the form of the solution in the vicinity of bifurcation points. 
From the equation for the kernel of the corresponding linear operator 




w"" + (p 2 - T)w" + Kw 



V 



u 
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and eq. (24) the relation determining position of bifurcation points is obtained: 

K l 1/2 

Pn = (nn) 



1 + 



(mr) 2 (mr 



nn + l-^ + l-^ + tiW, n = 0,1,2,.... (41) 
2 (mr) 2 {rm) A 



Clearly if the spring constant is small, then the nontrivial bifurcation 
points are located at same positions p n as in the Thurman and Mote model. 
The effect or elastic support is important for very small n (e.g. n = 1 or 
n = 2) and large values of K. However, if K is large, and (l/2)r/wr small in 
comparison with (mr), their locations may be distributed as in Fig. 5. A brief 
analysis of the above expression shows that the first nontrivial bifurcation 
point and some of the subsequent ones (usually for the first few n's ) are 
to a large degree determined by the value of K, and the same applies to the 
distances between the bifurcation points. Afterwards the term (mr) dominates 
and the bifurcation points are equally distributed. A complete nonlinear model 
consisting of equations (19) and (20) and including the term corresponding to 
the elastic support is very interesting both from mathematical and physical 
aspects; however, due to its complexity it will be analyzed elsewhere (Rajkovic, 
2005). 



3 The complete nonlinear model of Paidoussis 



3.1 Symmetry considerations 



The derivation of the complete nonlinear model, based on the work 
(Semler et al., 1994), in its correct form is given in Paidoussis (2003). This 
model, as presented in Section 2, in contrast to the previously analyzed mod- 
els, includes effects of gravity and pressure at the downstream end, so the 
dimensionless equations of motion of an extensible cylinder conveying fluid 
have the form: 



w"" + (p 2 - (r - u))w" + W + (r - a - n) (W + u 'w" + ^w' 2 w") 

-(3u"'w" + 4u"w"' + 2u'w"" + w'u"" + 2w" 3 + 2w' 2 w"" + 8w'w"w"') 



-7 



w'u' + ^w' 3 -{1-0 (-w" + u"w' + u'w" + \w l2 w"^ 



0, (42) 
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2 W 



.12 



i)w' w '' 



(p 2 - A)u" - {w"w'" + w'w"") - 7 

+ (r - a - n) w'w" = o. 

The relevant quantities have been defined in the previous section. 



(43) 



To begin with, it is of interest to inspect the symmetry properties of this 
equation. It is immediately clear that in the first equation, only odd powers 
of w turn up, while this is not the case with the second equation. Hence, 
solution of the set (42) and (43), v = (w u) T , satisfies the following symmetry 
condition 



(44) 



where T represents the operator of the symmetry group. In the vicinity of 
bifurcating solutions, the solutions of the complete model have the form 





f(0 
«(0 



= Mp) 




(45) 



where wo(x) represents the solution of the linear equation 

w"" + [(p 2 - (r - n)) - 7 (i - i)] w " + 1W > = o. 

Acting with the symmetry operator T on equation (45) one obtains 



(46) 



w(0 




Hence, the bifurcation equation satisfies the relationship 

g(-^,p) = -g(v,p) 

so that g possesses the Z 2 symmetry. The above symmetry properties indicate 
that the bifurcation is necessarily of the pitchfork type (Golubitsky, 1985). The 
physical representation of this symmetry is a reflection across the longitudinal 
pipe axis. The use of only one equation in the model of Holmes, the one 
involving w, may be justified based on these symmetry considerations. 
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3.2 Critical velocity 



With the fluid velocity as the bifurcation parameter, the procedure for 
obtaining critical velocity values corresponding to bifurcation points is almost 
the same as in the model of Holmes. Specifically, equation (46) 

w »» + [( p 2 _ r + n) - 7 (i - t)] w » + 1W > 
= w »» + (p 2 - r + nK - 7 ^[(i - 0w'\ = o. (47) 

corresponds to the linearized version of equation (8) of the Holmes' model. The 
only difference is the inclusion of the pressure term. The equation analogous 
to equation (B.8) in Appendix B determining bifurcation points is: 

2 (C 2 + 7) 3/2 - C 3 

= rar, 

3 7 

where 

C 2 = p 2 -(r-H)- 7 . (48) 
Taylor expanding and retaining terms to second order yields the following 
expression for the critical velocity values (bifurcation points): 



p n (L) =nn + -^- + -± 1 + 0(2) 

4 rnr 2 titt 

= n7r + ^ 7 + 2(r ~ n) +0(2), n = l,2,... (49) 
4 nir 

An important feature of the critical velocity p n , in contrast to previously 
considered less complex models, is its dependence on the pipe length through 
the length dependence of T and n. This dependence clearly diminishes with 
increasing velocity. Expression (49) is obtained assuming that fluid flow is in 
the direction of gravity (7 > 0). Bifurcations occur for velocity values higher 
than in the case when gravity, pressure and tension are not taken into account, 
under the assumption that 

7 + 2r > 2n. (50) 

The influence of gravity and tension is diminished by the pressure at the 
downstream end as intuitively expected. The individual effect of the gravity 
term becomes dominant if the pressure value approaches the value of tension. 
If the flow is in the direction opposite to the direction of gravity (7 < 0), the 
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expression for bifurcation points is 

ft>(L)=fW _ijL + i£zia +tf(2) 

Ann 2 mi 

i7-2(r-n) 

= nn-- + ^ - + ^2, n=l,2,... 51 

4 nn 

Clearly, gravity in this case acts in the same direction as pressure, and together 
they oppose the effects of tension. Additional insight into the position of bifur- 
cation points, in units of velocity, may be obtained by considering expressions 
for interpoint distances. For 7 > the interpoint distance is 

a ,rs 1 7 1 r-n 

Ap n (L) = n- 



Ann(n + 1) 27rn(n + l) 

= *-- A — r — 7r[7-2(r-n)]. (52) 
47rn(n + 1) 

As evident from the above expression, at low fluid velocities this distance is a 
sequence of increasing values provided that 

7 + 211 > 2r, 

and a sequence of decreasing values for 

7 + 2IT < 2r. 

Since 7 is length dependent [Eq. (2)], the former condition is more easily 
satisfied for long pipes, while the latter is more likely fulfilled for short pipes. 
For large n (high velocities), distance between bifurcation points has a fixed 
value of 71. For 7 < 0, expression (52) becomes 

. / r \ 1 7 1 r-n 

A Pn (L)=7T + 



47rn(n + l) 27rn(n+l) 

= + - — r — rrb - 2(r - n)]. 

4vrn(n + 1) 

In this case the distances between the bifurcation points form a decreasing 
sequence if 

7 + 211 > 2r, 
and an increasing sequence provided that 

7 + 211 < 2r. 

Recalling that 7 < 0, the first condition is more likely in short pipes while the 
second one is more probable in long pipes. 
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3.3 Normal form of the bifurcation equation 



The normal form of the pitchfork bifurcation modulo higher order terms 
(hot), reads 

9(P, p) = + 9wPP + hot. (53) 

In Appendix B, terms and g^ p are defined with the symbol x replacing 
li. In previously considered models the term g w was explicitly determined 
and in all cases it was negative. Hence, the sign of g mil determines whether 
bifurcation is of supercritical or subcritical type. For the complete nonlinear 
model, evaluation of terms g MWt and g w is much more complicated, and in 
order to minimize computational effort it is sufficient to determine just the 
sign of g^p. Following evaluation of the expression for g w given in the Appendix 
A, the following relationship is obtained: 

g w = (v*, L p ■ v ) = 2p (vq, Vq) , (54) 

where L is the linear operator corresponding to the set of equations (42) and 
(43) 

L = ( < + Kp 2 - (r - n)) - 7(1 - OK + Wo 
{ < 

Vq represents the solution of L = evaluated using procedure presented in 
Appendix B. With £ given by expression (48) v® has the following form 




v 



MO 
«o(0 



-J. 



(C 2 + 70 



'2^ 
,3 7 



J2/3 o~ J - 2 / 3 



'2(C 2 + 70 



3/2 ~ 



7 



2/3 



'2C 



n 3 7 J 2/3 \3 
satisfying boundary conditions 

v (0) = v (l) = 0. 
Explicit expression of the inner product in (54) yields 

5V = 2p ( v v'od£. 

J 



'2 (C 2 + 70 3/2 ^| 



7 



(56) 
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Performing integration by parts and using boundary conditions, one obtains 

g w = 2p[(v v' )\ 1 - j ' v%d£ = -2p J ' v%d£. 
Since the velocity p > and 

/\ 2 ^>o, 

J o 

it follows that 

9w = -2p fvldi<Q. (57) 

J 

Hence, once the sign of g^p is known an explicit evaluation of this term is 
unnecessary. However, an explicit evaluation of the term g MW is essential in 
order to obtain conditions determining the bifurcation type. 



Details related to the evaluation of g MMAt , due to its mathematical com- 
plexity and extensiveness, are presented in Appendix C. Only results relevant 
for the final form of the expression determining the sign of are presented 
in this section. Following appropriate calculations g^p assumes the following 
form: 



5Wm 



+A(C) 




(5* 



Terms A(£), A(£) and fi(£) are polynomial functions of ( and their explicit 
representation may be found in Appendix C. As in previous, less complicated 
models, two cases may be considered. 



Case 1: ^y/( 2 < 1. In the high velocity limit p 2 n — A > 0, so that the 
dominant part of (58) determining the sign of <7 MW1 is 

-2A(C)(r-.A-n + 7). (59) 

Hence, the bifurcation is supercritical or subcritical if the above expression 
is greater then or less then zero respectively. In order to be precise, it is 
informative to consider expression A(£). Explicitly, this expression is 
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where a and b are positive constants: 




and T (.) is the gamma function. The sign of A(() is also positive as 




irrespective of the possible negative signs of either J 2 /3 or J_2/3- Therefore, 
conditions for supercritical and subcritical bifurcations are determined by the 
term — (r — A — H + 7), so that conditions are: 

Supercritical condition: 

r + 7 <^ + n, (61) 



Subcritical condition: 

r + 7 >^ + n. (62) 

Recalling expressions for dimensionless system parameters (2), (3) and (21), 
inequalities (61) and (62) may be also expressed as: 

T(L) + (M + m)gL< EA + P(L), (63) 

and 

T(L) + (M + m)gL > EA + P{L). (64) 

revealing their dependence on the length of the pipe. Physically, the effects of 
gravity are related to the length of the pipe [equation (2)], so larger 7 may be 
associated with a longer pipe and smaller 7 with a shorter one. For short metal 
pipes 7 is rather small and its effects in inequalities (61) and (62) are weak. A 
comparison with corresponding inequalities for the Thurman and Mote model 
which includes curvature effects, shows that the pressure term replaces the 
velocity-dependent term. The reason is that the pressure term acts the same as 
the velocity term and it is clear that an adequately high level of pressurization 
alone may cause supercritical bifurcation. For pipes made of elastic material, 
gravity effects are important, as are the effects of axial flexibility, so that a 
supercritical bifurcation is more probable in shorter pipes. Intuitively, short 
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pipes acquire a buckled shape in an evolutionary manner (corresponding to the 
supercritical case), while longer pipes suddenly deform (corresponding to the 
subcritical case). For long pipes tensioning term is small (internal dissipation of 
the pipe material is assumed to be of the Kelvin- Voigt type) so that subcritical 
bifurcation is more likely to occur for low levels of pressurization. Thus, high 
pressurization in short pipes makes supercritical bifurcation more probable, 
while increasing the length of the pipe enhances the probability of subcritical 
bifurcation. 

Case 2: C 2 /7 < 1- Assuming that p 2 n — A < in this low velocity case 
(which is of real physical importance) the sign-dominant term is 



if + (r - a - n + 

where ip and tp are the following expressions: 

2^ 




(65) 



(66) 



and where Ci, c 2 and c 3 are constants (given in Appendix C). Clearly ip > 
and it is easy to notice that the sign of ip < does not depend on signs of 
J-2/3 ((2/3)C 3 /7^ and J2/3 ((2/3)£ 3 /"A , since the term containing these two 
Bessel functions is much smaller then the other two terms. Hence, the following 
conditions are obtained for bifurcation types: 

Supercritical condition: 



r + 7 >^ + n + ^, 



(67) 



Subcritical condition: 

r + 7 <^ + n + M. (68) 

Compared to conditions (61) and (62), inequalities pertaining to the low ve- 
locity case have an additional term acting in the same manner as the pressure 
and axial flexibility. The term \ip\ /ip is approximately proportional to 7 so it 
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counterbalances the gravity term on the left side of inequalities and the above 
inequalities may be further replaced with expressions 



T(L) > EA + P(L), 



and 



T{L) < EA + P(L), 



for subcritical and supercritical cases respectively. Since tension is inversely 
proportional to the length of the pipe, for long pipes subcritical condition is 
practically always satisfied. For short pipes supercritical bifurcation is possible 
only if pressurization term is small compared to axial flexibility, since the 
corresponding condition may be written as 



In the above expression a viscoelastic material has been considered, and since 
we are considering a time independent model, axial tension is 



The other possibility for this low velocity case, p\ — A < 0, has no physical 
meaning since n is small. 

Although a subtle interplay between system parameters requires care- 
ful analysis under specific circumstances of interest, in general it may be con- 
cluded that in the case of high fluid velocity a supercritical bifurcation is more 
likely in a short pipe, while a subcritical bifurcation is more probable in a long 
pipe (provided that appropriate inequalities (61) and (62) are fulfilled). In the 
case of low fluid velocities (along with (? n > A) the situation is the oppo- 
site: a supercritical bifurcation is more likely in a long pipe, while subcritical 
bifurcation is preferred in short pipes provided condition (69 ) is satisfied. 



4 Conclusion 

Exact analytical solutions in the vicinity of bifurcation point for each 
model are obtained along with the derivation of conditions that classify bi- 
furcations as supercritical or subcritical. Moreover, the analysis is performed 
in such a manner that influence of important quantities on the dynamics of 




(69) 



T = a A = (Ee + E*e)A = eEA. 



where averaged axial strain e due to lateral deflections w is 
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supported fluid conveying pipes may be analyzed in the light of increasing 
complexity of each model. 

Two important features of the stationary bifurcations for the supported 
fluid-conveying pipes should be emphasized. First, all bifurcations are of the 
pitchfork type as a consequence of reflection symmetry. Second, all pertur- 
bations of the pitchfork bifurcation, due to gravity or curvature effects for 
example, preserve the topological form of the unperturbed bifurcating dia- 
gram, due to the fact that remains the solution of the perturbed equation. 
Consequently, unfolding of the bifurcation does not take place which requires 
that is not the solution of the perturbed equation. 

Models that consider both axial and lateral deflections, hence two- 
equation models, enable the possibility of both supercritical and subcritical 
pitchfork bifurcations. In contrast, the single equation model of Holmes which 
considers just transverse deflections allows only supercritical bifurcations. An 
important general characteristic of classification of generic codimension-1 bi- 
furcations into supercritical or subcritical is that the bifurcation type depends 
on only two factors: nonlinear terms and boundary conditions. Hence, non- 
linear terms figuring in the equation for axial deflections make an important 
contribution to the terms defining the normal form of the bifurcation equation. 

In the complete nonlinear model of Paidoussis critical velocity values 
at which bifurcation occurs depend on tension, gravity and pressure. For 7 > 
gravity and tensioning exert matching effects reflected in shifting bifurcation 
values in the positive direction with the shift due to gravity being one-half 
the corresponding tensioning shift. The effects of pressurization oppose effects 
of gravity and tension. For 7 < 0, gravity and pressurization act in the same 
direction, while tensioning exerts opposing effect. The critical velocity values 
for each model are presented in Table 1. A summary of conditions classifying 
bifurcations as supercritical and subcritical is presented in Tables 2 and 3. 
An increase of complexity of the bifurcation type conditions may be traced, 
starting with the least complex model of Thurman and Mote and ending with 
the complete nonlinear model of Paidoussis. The velocity-dependent term ap- 
pearing in inequalities corresponding to the Thurman and Mote model with 
curvature transforms into an analogous pressurizing effect in the most complex 
model. Setting 7 = n = 0, and assuming that T and n are independent of the 
pipe length in the complete nonlinear model of Paidoussis, one obtains classi- 
fication conditions for the Thurman and Mote model. Although conditions for 
determining whether bifurcation is of supercritical or subcritical type involve 
a delicate interaction among system parameters (gravity effects, tensioning, 
pressurization and axial flexibility), it may be concluded that a general ten- 
dency is that for high fluid velocities supercritical bifurcation is more likely 
in short pipes, while subcritical is more likely in longer pipes (provided that 
appropriate inequalities (61) and (62) are satisfied). When the fluid velocity 
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is low, but still high enough that the square of it exceeds axial flexibility, a 
supercritical bifurcation is more likely in longer pipes (with conditions (67) 
and (68) in effect), while a subcritical bifurcation may occur in short pipes 
provided that condition (69) is satisfied. 



5 Acknowledgement 

The authors wish to thank anonymous referees for helpful and instruc- 
tive comments, suggestions and criticisms. This work is partially supported 
by the Serbian Ministry of Science and Technology as part of the project 01 
1986. 



24 



Appendix A. Lyapunov-Schmidt reduction 



Reduction of a nonlinear equation or a system of nonlinear equations 

$ t (y,A) = 0, i = l,...,n (A.l) 

to the single algebraic equation g(x, A) is the essential feature of the Lyapunov- 
Schmidt procedure. The vector y = y(yi, ■■■,y n ) is the unknown in the above 
equation, while A is a vector of parameters. We assume that only one param- 
eter A is of concern. The main starting assumption is 

$,(0,0) =0, (A.2) 

and of interest is to describe the solutions of this system locally, in the vicinity 
of the origin. If the rank of the n x n Jacobian matrix L = (rf$)o,o, is equal to n 
(nondegenerate case), the implicit function theorem guarantees the existence 
of solution y as a function of A. If the rank is not equal to the size of the 
Jacobian matrix, we assume the minimally degenerate case, i.e. 

rank L — n — 1. 

Assuming that $ : R n x R k -> W 1 is a smooth mapping, vector space comple- 
ments M and N are chosen to ker L and range L, respectively, so that 

R n = ker L © M, 

and 

R n = N © range L, 

where denotes the null space of L. Introducing the projection operator 
E : W n — > range L, the starting system of equation is expanded into 

£$(y,A)=0, (A.3) 
(/-£)<%, A) = 0, 

where I— E is the complementary projection operator to E. Solving the first 
equation of (A.3) which fulfills the conditions of the implicit function theorem 
for n — 1 of the y variables, and inserting solutions in the second equation, 
yields an equation for the remaining one variable. Because of the splitting of 
]R n , any vector y G K" may written as y = v(e ker L) + w(e M), so that 
mapping 

F : (ker L) x M x R k -> range L, 
is given by expression 

F(v,w, A) = E$(v + w,\). 



25 



The linear map 

L : M — > range L 

is invertible, thus according to the implicit function theorem it may be solved 
for w near the origin. Denoting this solution as w = W(v, A) : ker Lx IR fe — > M, 
which satisfies 

E®(v + W(v,X),X) = 0, W(0,0) = 0, 

and it may be substituted into the second equation of (A. 3) to obtain the 
reduced mapping : ker LxK'^ N, where 

<f)(v, A) = (/ - E)$(v + W(v, A), A). 

Consequently, the zeros of the <j)(v, A) are in one-to-one correspondence with 
the zeros of A), or explicitly 

<j)(v, A) = if and only if $(u + W(v, A), A) = 0. 

The reduced function (f)(v, A) may be further used to obtain the algebraic 
equation 

g(x,X) = (vq, (j)(xv , A)) , 

where Vq G (range L)^, and (., .) represents the standard inner product. Illus- 
tration of the reduction procedure is shown in Fig. Al. Finally, the derivatives 
of g figuring in the Taylor expansion of g(x, A) around the origin, after com- 
putation in terms of the original mapping A), are given below: 

9x = 0, 

9xx = (vQ,d 2 ^( y v ,v )) , 

9xxx = (vq, d 3 ^(v , v , v ) - 3d 2 $(t; , L~ 1 Ed 2 §(v , v ))) , 

gx = (v* ,$x), 

g\x = (Vq, d$\ ■ Vq - d 2 $(vo, L^EQx) ■ 



Appendix B. Gravity effects in the model of Holmes 

Solution of equation (11) 

t V + tw' + (t 2 - (\) 2 )w = 0, (B.l) 

may be expressed as 

w(t) = CiJi/ 3 (t) + C 2 J-l/ 3 (t) 
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so that solution of eq.(ll) may be written as 



u{t) = 



7 



1/3 



ClJl/3 



'2a 2 + 7 3/2< 



+ c 2 J_i/3 



, 3 7 
3 7 , 



(B.2) 



Bessel functions posses the following well known properties, cf. Gradshteyn 
and Ryzhik (1994): 



d ( J u (x) 
dx 
d_ 
dx 



X" 



so that 



x 



U+1 J u { X )dx = -X-^Jv-^x), 

J x v+1 J u (x)dx = x u+1 3 u ^(x). 
Furthermore, boundary conditions of Eq.(ll) 



u 



k=o= 



yield 



u' k=i= 0, 



/2« 3 \ /2 a 3 

ClJ-2/3 [ 3— 1 -C 2 J 2/3 (- — 



0, 



ClJ-2/3 



'2(a 2 + 7 ) 3/2 ^ 



7 



+ c 2 J 2 /3 



3 7 

2 ( a 2 +7 )3/2- 

3 7 , 



0. 



The condition for obtaining nontrivial solutions c±, c 2 7^ requires 



J 2/ 3 fl^ J 



-2/3 



'2(a s 



\3/2 l 



J 



2/3 



,3 7 
'2 (a 2 + 7 ) 3/2 \ 

,3 7 ; 



J. 



2/3 



7 

/2a 3 
I 3 7 



+ 



0. 



(B.3) 



(B.4) 



(B.5) 



(B.6) 



Assuming weak gravitational influence (7 ~ 0), the arguments of the Bessel 
functions in the above equation are large so that asymptotic expressions for 
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Bessel functions may be used 



Hence, 



^(x) = J— cos(x - % - j) + $(x- 3 / 2 ), 

V 71X I 4 

J-„(a;) = \ — cos(x + \v - + tf(ar 3/2 ). 
V ttx 2 4 



' 2/:i( ' r) = \l^x~° 0S ( X ~ T^' 



/ 2 77T 

J-2/3(x) = J— cos(x+— ). (B.7) 



In a straightforward manner, condition (??) yields 



2 (a 2 + 7 ) 3 / 2 



3 



nvr. (B.8) 



3 7 

Recalling that a 2 + 7 = p 2 , the above condition is equivalent to 

P 3 " (P 2 - if 2 = t^T, 

which, assuming 7 << p may be Taylor expanded, and retaining terms to 
second order yields 

Keeping in mind the substitutions introduced in order to cast the equation in 
an analytically solvable form, the solution we seek is 

H0= with«(0) = 0. (B.9) 

Jo 

Inserting eq.(70) in (B.3) the following expression is obtained: 

= ^ (^) 2/3 i% 2/3 J./ 3 (»)* ( + $ (^) 2/S f>J-l/ 3 (»)*7, 

(B.10) 

so that 



+C 2 



(a 2 + 7 0J- 2 /3 (| 



2(a 2 + 7 3/2 ' 



a 2 J. 



7 



•2/3 



(a 2 + 7OJ2/3 



'2 (a 2 + 7O 



3/2 ~ 



7 



-a 2 J 



2/3 



/2a^ 
V3 7, 



(B.11) 
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Constants in the above two expressions satisfy the following relations: 



c 2 _ J-2/3 (l^rj c 2 



,3 7 



Cl 



C l J 2 /3 (I 

Finally, the complete solution v(x) may be written as 



(B.12) 



h/s {3- J - 2/3 U — 



2 a 3 



- J - 2/3 I 37 J J2/3 I 3 



2 (a 2 + 7 



3/2^ 



7 



7 



for 7 ~ 0, this acquires the asymptotic form 

1 



«(0 - 



(« 2 + 70 : 



sin 



2 (a 2 + 70 

3 7 



3 / 2 - a 3 



while 

as expected. 



limw(0 = sin(n7r0, 

7^0 



(B.13) 



(B.14) 
(B.15) 



Appendix C. Normal form of the bifurcation equation for the com- 
plete nonlinear model 



The normal form of the pitchfork bifurcation modulo higher order terms 
(h.o.t), reads 

g(v> p) = ? w i" 3 + QwW + h -°- t - 

As shown in the main part of the paper, explicit determination of the term g Mp 
is not necessary since only its sign is relevant, and it was demonstrated that it 
is negative. The remaining term, g^^, requires explicit determination in order 
not only to evaluate its sign which enables classification of bifurcation into 
supercritical or subcritical type, but also to extract conditions, in the form 
of inequalities, that need to be fulfilled in order for each bifurcation type to 
arise. In Appendix A, it was shown that this term requires evaluation of the 
inner product which may be written as 

9xxx = (v* \d 3 <S>{v ,v ,v )) - (u5|3^$(uo,£~W$(^,uo))- (Cl) 
Noting that the adjoint operator L* of 

T ( <" + [(p 2 - (r - n)) - 7(1 - OK + 1< \ 
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is equal to L, the following expression is obtained for the first inner product 
in (C.l): 

(v* \d 3 <S>(v ,v ,v )) = 9a J Q v' 2 v' \d£ + 97 y (1 - £)vgvZv d£- 

37 C v' 3 v dC + 12 f 1 v^v d^ + 12 v*v™v d£ + 48 /' > 
j j j j 

where 

a = r — .4. — n, 

and where v represents the solution of L = evaluated using procedure 
presented in Appendix B. With ( given by expression (48) v has the following 
form 




(C.2) 



Once the inverse operator L~ 1 Ed 2 &(v , v ) is evaluated following a lengthy 
procedure, the second inner product requires evaluation of the following inte- 
grals: 



'vl\?>d 2 <5>(v ,L- l Ed 2 <5>(v ,v {i )) = 

(<*o + 7) / W + v")vov'odt, - / £r)'voVod£- 
J J 

I" Zv"v v' dS - 3 P V '"v vZ - 4 f 1 rf'vov'''^ 
Jo Jo Jo 

-2 l\""v' Q di- 1 [\'v v' dC, 

Jo Jo Jo 



(C.3) 



where 



2(a + 7) 
p 2 -*4 



.A 



I I III! 1 // lll\ Jt- 



(C.4) 



Evaluation of integrals was performed using approximate expressions 
(??). Results were checked using lower and upper bounds for the Bessel func- 
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tions, given by the following inequality (Neuman, 2004): 



1 



x 



< 



X 



r(a + l) \2J 3(a+ 1) 



r(a + l) V2 
1 



x 



cos 



2(a + 1, 
2a + 1 + (a + 2) cos 



2a + 2 



(C.5) 



After lengthy calculations the following expression is obtained 



where 



HO 



(r - .a - n - 



'2£ 

3 7 



A(C) 



P^-^ V37/ 



2 Nl/3 

' J C + n(C) 



3 

2 \ 8 / 3 



A(C) = 

a and b are positive constants 



M^) ( j -/ 3 (^)) J 2 /s0^ 



-2/3 



'2C. 
3 7 



3\ \ 4 ' 



(C.6) 



(C.7) 




Furthermore, term A(£) is 



where 



a(o-*(I)V(^(I?))^(I?) + 



r(l)lr(i) r 
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Finally, term Q(() is equal to 

-(</? + (r- .4- n + 7 )v), 

where 



* = (r) 7 ( J - 



2/3 



3\ \ 2 



2C 



3\ \ 2 



C2 U 7 



2 \ 8 / 3 / ^2C 

J -2/3 



37 

3\ \ 4 



37 



3 7 

J-2/3 



3\ \ 3 



- c 3 ( -7 



-Ml?)) 4 



3 7 

2 \ 7 / 3 



'2/3 



'2£ 
3 7 



777 + o7 



2 i/3 



Ci, c 2 and c 3 are the following constants: 



ci = 12 



2 i/s 



r(J)r(J) Vr(i) r(|) 



C 2 



Vr(|) r(|)r 



2 l/3 



C 3 = 



Vr(|) r(i 
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Figure captions 

Fig. 1. A pipe with supported ends conveying fluid. The diagram also 
shows coordinates used in the text; p represents the fluid velocity. 

Fig. 2. Supercritical pitchfork bifurcation for the model of Holmes with 
and without gravity and tensioning effects. Fluid flow is in the same direction 
as gravity. 

Fig. 3. Supercritical pitchfork bifurcation for the model of Holmes. 
Fluid flow is in the direction opposite to the direction of gravity. 

Fig. 4. Subcritical pitchfork bifurcation for the model of Thurman and 

Mote. 

Fig. 5. Supercritical bifurcation for the fluid conveying pipe with and 
without elastic support. Dashed lines correspond to the case without elastic 
support. The case with elastic support assumes if»0, and (l/2)r/rar <C nir. 
The first three modes are shown. 

Fig. Al. A geometrical interpretation of the Lyapunov-Schmidt reduc- 
tion. 



35 



Tables 

Table 1. 



Model 


Critical velocity p n 




7 > 


7 < 


Holmes 


4 n-K 2 n-K 


4 ri7T 2 J17T 


TM and TM-k 2 (7 = 0) 


2 n-K 




TM-K (7 = 0) 


n7r +2(n7r) 1 2(n7r)3 




Complete nonlinear model 


4 n7r 2 mr 2 ri7r 


4 mr 2tik 2 nK 



Table 1. Critical velocity values for models of Holmes, Thurman and Mote 
(TM), Thurman and Mote with curvature (TM-K 2 ) and the complete nonlinear 
model of Paidoussis. 
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Table 2. 





High velocity limit 


Model 


Supercritical condition 


Subcritical condition. 


TM 


r <a 


r >a 


TM-k 2 


r <^+f (nvr) 2 


r >^+f (nvr) 2 


Complete nonlinear model 


r + 7 <^+n 


r + 7 >^+n 



Table 2. Conditions for development of supercritical an subcritical bifurca- 
tions in the high velocity limit 
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Table 3. 





LOW VELOCITY LIMIT 


Model 


Supercritical condition 


Subcritical condition. 


TM 


T <A-\(nnf 


T >A-\(nnf 


TM-k 2 


T <A+l(nn) 2 


T >A+l(nn) 2 


Complete nonlinear model 


r + 7 >A+u+ l -f 


r + 7 <A+u+ l -f. 



Table 3. Conditions for development of supercritical and subcritical bifur- 
cations in different models of fluid conveying pipes. Term \ip\ /ip, proportional to 7, 
is denned in eq.(66). 
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